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We propose an easy to use model to solve for interacting atoms in an optical lattice. This model 
allows for the whole range of weakly to strongly interacting atoms, and it includes the coupling 
between relative and center-of-mass motion via anharmonic lattice terms. We apply this model to a 
high-precision spin dynamics experiment, and we discuss the corrections due to atomic interactions 
and the anharmonic coupling. Under suitable experimental conditions, energy can be transferred 
between the relative and center-of-mass motion, and this allows for creation of Feshbach molecules 
in excited lattice bands. 



PACS numbers: 37.10.Jk, 34.50.-s, 03.75.Lm, 34.50. Rk 



Optical lattices form a suitable environment for high 
precision experiments with interacting atoms. Two 
atoms can be isolated from other atoms by placing them 
on a single lattice site, and many sites can be filled si- 
multaneously. While the lattice parameters such as depth 
and geometry can be adjusted via the laser field, the in- 
teractions can be tuned using a Feshbach resonance by 
applying an external magnetic field. Feshbach molecules 
in a lattice, created by sweeping the magnetic field over 
resonance, can be transferred into deeper vibrational 
bound states, for instance, by applying stimulated Ra- 
man adiabatic passage 

Precise values for relative interaction strengths of ru- 
bidium atoms were derived by studying coherent colli- 
sional spin dynamics in an optical lattice 0, Q ■ These 
high precision measurements provide challenges for the- 
oretical coupled-channels models based on current state- 
of-the-art interaction potentials [4[ . One may wonder for 
instance, at what level of precision it is possible to calcu- 
late interaction properties before the Born-Oppenheimer 
approximation, which is the underlying foundation for 
the potentials, breaks down. 

Before conclusions can be drawn on the limitations of 
two-body interaction models, one has to make sure that 
the correct comparison is made between theoretical quan- 
tities of such a calculation, and the measurements that 
depend on these interactions. For instance, one cannot 
always put the resulting scattering lengths of a two-body 
collision, defined in the limit of zero collision energy, as 
the on-site interaction in a Hubbard model. The diver- 
gence of the scattering length on resonance will give rise 
to physically unrealistic large energy shifts. A resonant 
interaction takes the two-body interaction in the uni- 
tarity limit, where the solution of scattering wave-func- 
tions are shifted over it/2 compared to non- interacting 
atoms, and one would rather use expressions based on 
the (energy) dependent scattering phase shift. This argu- 
ment also applies for high-precision experiments on non- 
resonant systems, since small energy- dependent correc- 
tions already can be of importance. Also, the relative and 
center-of-mass motion of two interacting atoms, which 
can become coupled due to different atomic species [!, 0| 
and anharmonic terms in the lattice potential Q, can 



give rise to shifts in the on-site interaction. 

The model we put forward in this paper is conception- 
ally simple and easy to use. It is based on first-order 
perturbation theory starting from an existing solution of 
two interacting atoms in a harmonic potential. We show 
that the model is valid for moderately deep lattices, i.e. 
where tunneling to neighbouring sites is negligble. In this 
way, we are able to make a proper comparison between 
the high-precision measurements by Widera et al. @ and 
accurate rubidium scattering models. We demonstrate 
the importance of energy-dependence in the scattering 
phase shift and of anharmonic corrections for this exper- 
iment, and also show how experiments using a Feshbach 
resonance can produce molecules of a mixed relative and 
center-of-mass motion nature. 

This paper is outlined as follows: First we give a de- 
scription of our model. Then we apply it to a spin- 
dynamics experiment in a lattice. In the third section, we 
discuss the nature of the molecules that have a mixed rel- 
ative and center-of-mass motion, predicted by our model. 
We end with a concluding section. 



I. MODEL 

In the ultracold regime only s-wave interactions are 
allowed, and we conveniently model the interaction with 
a pseudo potential Q 



V int (r) 



aS (i) (r) — r. 



in 



dr 



(1) 



Here r is the distance between two colliding atoms, and 
a the s-wavc scattering length @. The use of this zero- 
range pseudo potential is allowed, since the range of the 
real interaction potentials is much smaller than the lat- 
tice spacing. The total Hamiltonian for a pair of inter- 
acting atoms in an optical lattice is the described by 



H = — 



1. 

2m V 2 



+ 2V; 



(2) 



+Mat (r+|) + VLt (R - I) + V hlt (r). 
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Here 

3 

MatW =V ^sin 2 (fc L a; J ) (3) 
i=i 

is the optical lattice potential with depth Vb and spac- 
ing d = 7r/fcL, where is &l the wavenumber corresponding 
to the laser frequency and x = (x%, X2, xs) is decomposed 
along the lattice axes. R and r denote the center-of-mass 
(CM) and relative position of the atoms. The restriction 
to a pair of interacting atoms is an approximation, but 
exact for deep lattices (for example in the Mott-insulater 
regime) for two trapped atoms per lattice site. Separa- 
tion of the CM and relative motion is possible when each 
site is treated as an harmonic oscillator (HO) with fre- 
quency lu = \/2V{)kf j /m. The lattice potential for two 
particles can then be written as 

Mat (R + F -) + Mat (R - |) « im^ 2 (2R 2 + l -v 2 ). (4) 

An exact solution for two interacting atoms in a HO trap 
can be derived for the relative motion, as was first shown 
by Busch et. al [§]. However, lattice effects beyond the 
harmonic approximation, i.e. anharmonic terms as well 
as tunneling to neighbouring sites, can be of significant 
importance even when the atoms are considered trapped 
in a single site ■ The consequence of the introduction 
of such terms is that the CM and relative motion become 
coupled. To tackle this problem, we handle anharmonic 
and tunneling terms in a pcrturbativc procedure, starting 
from the separable solution: 

4° ) (R,r) = $ s (R)^(r). (5) 

Here $g(R) and (j) s (r) are the exact solutions of the CM 
and relative motion, corresponding to states labeled with 
S = {Si, S2, S3} and s = {n, l,m}, respectively. We 
use different quantum numbers for CM and relative mo- 
tion, to adapt optimally to both the symmetry of the 
perturbation term and the symmetry of the interaction 
region. Therefore, CM is always described in cartesian 
coordinates, with quantum numbers Sj = 0,1,2,.. for 
the different lattice axes. Relative motion is decomposed 
in spherical coordinates, with principal quantum num- 
ber n = 1,2,3,.., and orbital quantum numbers / and 
m. Without the presence of atomic interactions, i.e. 
when a = 0, many HO states are degenerate. The corre- 
sponding collection of quantum numbers is denoted with 
the symbol {■, •, ■}. 

For the case a — 0, one can also find an exact lattice 
solution completely in cartesian coordinates, for example 
in terms of Mathieu functions on single-particle coordi- 
nates. This allows us to estimate the importance of tun- 
neling compared to anharmonic terms, shown in Fig. [1] 
Tunneling gives rise to a broadening W — E t — E b of 
the HO level structure, which results in the band struc- 
ture of a square lattice [Tcj, and also has a level shift 
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FIG. 1: Estimation of tunneling and anharmonic effects based 
on exact lattice solutions for non-interacting particles. The 
ratio of band width W and level shift AE is shown as 
function of lattice depth Vo, given in units of recoil energy 
E TCC — h 2 kf j /2m. The black solid and red dash-dotted lines 
correspond to the lowest and second but lowest symmetric 
bands, respectively. The blue dashed line represents the low- 
est anti-symmetric band. 

AE = Eno - (E t + E b ) /2. Here E t and E b are the band 
top and band bottom energies, respectively, and -Eho is 
the corresponding HO level. 

We find that sufficiently deep in the Mott phase, in 
particular for lattice depths as used in current experi- 
ments, that anharmonic terms dominate above tunnel- 
ing: W < |AE|. Since also \AE\ < E w , we expect 
in the case of non-zero interatomic interactions, a per- 
turbative approach with anharmonic corrections only is 
sufficient to yield accurate results. In the remainder, we 
therefore restrict ourselves to single-site solutions includ- 
ing anharmonic terms. As perturbation term we then 
have: 

H'(R, r) = Vi at (R + + Mat (R-0 + (6) 

-imw 2 (2R 2 + ir 2 ). 

We define our perturbative solutions in accordance with 
standard perturbation theory: 

s Ss = 4 s ) + <4 iii/'i4°i) + ..., (7) 

and 

, _ . (o) , v (Vflavg) .(o) , f8 , 

VSs - IPss + 2^ p (0) _ p (0) *W + ■ ■ ■ ' ^ 
s'fts ^Ss ^S's' 

for the first order corrected energies and wavefunctions, 
respectively. E^ = Es+e s is the sum of the HO energies 
corresponding to the CM and relative motion. To check 
the accuracy of our model, we compare the results with 
the exact lattice solution available for the case without 
interaction. Regarding the energy, we compare E$ s with 
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FIG. 2: Comparison of the exact and perturbative wavefunc- 
tions in the lowest band, for a = 0. Shown are the zeroth 
order HO solution (dashed line) , the first order corrected solu- 
tion (red solid line) and the exact localized Wannier function 
(dash-dotted line), for CM coordinate R = 0. The relative 
coordinate is scaled on aho, with a^ — ^/2fi/(mw) the har- 
monic oscillator length of relative motion. 



-E-cxact = {Et+Ejj) /2. With only first order corrections we 
already have (Es s -E eKact )/E exac t < 0.7%, for the lowest 
band at a lattice depth Vq = 25E tec . Note that the accu- 
racy is of the same order as the relative contribution of 
tunneling (Fig. [I]), and therefore the meaning of second 
order terms is limited. However, this first-order accuracy 
is sufficient for our description of the spindynamics ex- 
periment in Section III. Naturally, the accuracy improves 
when lattice depth is increased. We also compared the 
results for the wavefunction, and an example is given in 
Fig. O for the same value of the lattice depth. Shown 
are the first order corrected wavefunction, together with 
the zeroth order HO solution and the Wannier function, 
being the exact solution for a single site. Anharmonic 
corrections make the trap less tight compared to the HO 
approximation, resulting in a small decrease of the prob- 
ability for finding the particles near r — around the 
origin, corresponding to a small increase of probability 
density in the barrier. The difference between the Wan- 
nier function and the HO wavefunction with anharmonic 
corrections is small compared to the difference between 
the latter and the unperturbed HO wavefunction. This 
shows that the method converges quickly. 

When we include interactions (a^O), the zeroth order 
relative wavefunction and energy depend parametrically 
on the scattering length according to the Busch model Q : 
4> s (r) = <f> s (r;a) and e s = e s (a). Hence the anharmonic 
correction becomes a function of scattering length as well. 
Note however that the perturbation term itself does not 
depend on the scattering length, illustrating that the in- 
teraction is modeled exactly. 

Finally, we show how to construct a solution for the 
whole lattice. Due to the relative interactions, there is 
no periodicity for the relative motion. However, the CM 
motion remains periodic. In fact the CM boundary con- 



ditions are dictated by the Bloch-theorem. These bound- 
ary conditions are satisfied when we construct the lattice 
solution as 

^ Q 5 S (R,r)=^e i Q- D ^ Ss (R-D,r), (9) 

D 

with direct lattice vector D = (tt-ic?, n^d, n^d), integers 
rii, and Q the CM quasi momentum. This demonstrates 
that the wavefunctions ipSs obtained with our model 
are effectively Wannier functions of a pair of interacting 
atoms in an optical lattice. 

II. SPINDYNAMICS 

In this section we apply our model to the spindynamics 
experiment, as carried out in the Bloch group [2, l3j. This 
allows us to investigate systematically the effects of an- 
harmonic terms and a non-zero interatomic interaction. 

In an optical lattice one is able to trap several spin- 
states at the same time. Starting with atoms prepared 
in a specific one-particle (hypcrfmc) spin-state |/,m/), 
collisions between two such atoms gives access to other 
two-particle spin configurations. When only weak mag- 
netic fields are applied, the total magnetization is con- 
served, and therefore coherent collisions between two- 
particle states of equal total two-particle spin F occur. 
This can be described by a Rabi-like model. For atoms 
prepared in / = 1, which is the case we will treat here, 
this is only a two level system, with effective Rabi fre- 
quency 

njf = y/s^ + sff. (10) 

Here fijf is the bare Rabi frequency depending on the 
coupling strength of the spin-changing collision. Detun- 
ing <5if contains two contributions: 

Si^So + SiB 2 ). (11) 

So is given by the difference of two interaction energies 
corresponding to collisions that leave the spin configura- 
tion unchanged, and 8{B 2 ) is a second order Zeeman shift 
between the initial and final states. By performing the 
experiment at different magnetic field strengths B, the 
Bloch group was able to substract flu and 5q, thereby 
being able to derive precise values for relative interaction 
strengths of rubidium. Treating the interaction energy 
as a linear perturbation in the parameter k^a, one can 
derive the Rabi parameters as 

2 v^2 1 
fiif = -^rU(a 2 - ao), 8 Q = —U(a 2 - a )- (12) 

The factor U is here defined as 

U= — x [ d 3 x|V| 4 , (13) 
m J 
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with ip the lowest HO eigcnfunction. U depends on 
the lattice depth, but is independent of the scattering 
length. Hence in this approach, differences between in- 
teraction energies are due to scattering-length differences 
only. The values of the scattering lengths ap, F = 0, 2, 
corresponding to collisions in subchannel F, are calcu- 
lated based on highly accurate rubidium potentials [J], 
and are given in Table |U Note that V 

can deviate signif- 
icantly from the above-mentioned solution of two inter- 
acting atoms in a trap, and therefore it can be expected 
that a proper treatment of interatomic interactions could 
have a large impact on the Rabi frequency. While this ef- 
fect will be mostly pronounced when close to a Feshbach 
resonance, an effect can also be expected when highly 
accurate measurements are performed, as in this experi- 
ment of Widera et al. 

The experimental results for the scattering lengths, 
based on the above Rabi model Q, agreed just within 
error bars with the predictions. This apparent discrep- 
ancy was most clearly seen for the / = 1 case. We will 
now investigate this / = 1 experiment by calculating the 
interaction strengths from our lattice model, and analyze 
the effects of anharmonic terms and exact interatomic in- 
teractions. Note that anharmonic corrections were also 
taken into account in Ref. Q , which already led to a bet- 
ter agreement between theory and experiment. We can 
most clearly compare theory and experiment by regard- 
ing the frequency at B = 0, since the field depen- 
dence does not depend on the interatomic interactions. 

As a starting point for comparison, we apply the Rabi 
model as described above, with U calculated from HO 
solutions, and we use the scattering lengths from Ta- 
ble U to obtain the effective Rabi frequency. This gives 
0|[ IO (B = 0) = 2tt x 49.27Hz. Note that for a com- 
parison between theory and experiment, it is sufficient 
to consider Q'^ at zero magnetic field, since the mag- 
netic field dependence does not involve any knowledge 
of the interatomic interactions. Then, we first calculate 
f2| f (-B = 0) by using Wannier functions in the expres- 
sion for U, in order to analyze the effects of anharmonic 
terms only. Wannier functions are exact solutions for a 
lattice without interactions. This is similar to the pro- 
cedure performed by Widera et al. [|[. Second, we want 
to analyse the effect of having exact interatomic interac- 
tions only, and calculate the Rabi parameters by using 
the solutions for two interacting atoms in a harmonic 



TABLE I: Theoretical predictions for the scattering lengths of 
the F — and F — 2 channel, for atoms with one-particle spin 
/ = 1, based on accurate rubidium interaction potentials Q. 
Here F is the total two-particle spin. The values are given in 
units of the Bohr radius ao. 



trap: 



where 



2\/2 



E iat (a)=E(a)-E(a = 0). 



(14) 
(15) 

(16) 



Here we define the total energy E(a) = Eh*] (a) — 
E s + e s (a), with S = {0, 0, 0} and s = {3, 0, 0}, accord- 
ing to the solution of two interacting atoms in a trap [ill ] . 
Third, the same is done, but also with the anharmonic 
terms included by taking E(a) = Es s {a), in order to com- 
pute the combined effect. The results are shown in Ta- 
ble[TTl by calculating the ratio of the different f2| f (B = 0) 
over the initial frequency fl'^°(B = 0). The table also 
shows a comparison with the experimentally obtained ef- 
fective Rabi frequency. In all calculations a lattice depth 
of Vo — 45-E' r ec is used. Note that the third calculation is 
the most precise one, with exact interactions, giving rise 
to modified wavefunctions and energy levels compared to 
the HO calculation. Moreover, the anharmonic effects 
are also included via perturbation theory up to a high 
precision, which follows from the table (discussion could 
also go to the end). 

From the results we find that anharmonic terms induce 
a negative shift of order 10 %. To the contrary, higher 
order interaction effects induce a positive shift, which is 
of order 3%. Hence, anharmonic corrections are domi- 
nant, while the first order approximation for the interac- 
tion energy is already fairly accurate. The net result is a 
7% improvement with respect to the initial model of the 
Bloch group. Although the theoretical and experimental 
values still differ by 30%, this is just within theoretical 
and experimental error bars. Note that the largest con- 
tribution to the theoretical error bar is due to the small 
difference 02 — ao, which is only a percent of the values 



TABLE II: (left column) Correction factor Q' i{ /Q'^° of the 
effective Rabi frequency compared with the same quan- 
tity computed in the HO approximation with the interaction 
treated as linear perturbation. The first row shows the re- 
sult when only anharmonic effects are taken into account. 
The result in the second row is obtained with only higher or- 
der interaction effects taken into account. In the third row 
the results are given of both effects acting together. Right 
column: corresponding values for the ratio between the the- 
oretical Q' if and the experimentally obtained effective Rabi 
frequency fijj Xp = 2tt x (35.4 ± 0.7)Hz. 



fflF [ap] 



101.78 ±0.2 
100.40 ± 0.1 



included corrections 




anharmonic 
interatomic interactions 
anharmonic + interatomic interactions 


0.897 1.25 
1.033 1.43 
0.935 1.30 
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of ao, a,2 themselves. A second issue is a possible system- 
atic error in U, related to the uncertainty in the lattice 
depth. At present there is no direct measurement of this 
coupling constant. 



III. ADIABATIC CREATION OF FESHBACH 
MOLECULES WITH A CM MOTION 

In the previous section we applied our model to non- 
resonant scattering of atoms in the lowest lattice band. 
This means that the energies of the relative and CM mo- 
tion are relatively close to their respective HO ground 
states. The coupling between CM and relative motion is 
taken into account, which also includes coupling to higher 
bands when second order lattice perturbation terms are 
incorporated. However, the interaction induced coupling 
between the bands is negligible since the energy distance 
to non-ground state CM motion is too large. When the 
interaction strength is increased, this is no longer true, 
and several relative and CM states become (nearly) de- 
generate, and therefore first order lattice perturbation 
terms will already result in an efficient coupling of these 
states. In this section, we show that this implies that 
atoms can be transferred into molecules, with a simulta- 
neous transfer of quantized energy of the relative motion 
to the CM motion, and vise versa. 

Tuning the scattering length through resonance can 
give rise to a transfer of atoms from one HO level to the 
next one Q, and under adiabatic conditions this can be 
observed experimentally [T^ |. In a pure harmonic trap, 
one is therefore able to increase the relative energy by 
2Hu>. This is illustrated in Fig. [3] by the thin gray lines. 
The system can become degenerate to other states which, 
in turn, can be lifted by anharmonic terms. In general, 
coupling strengths depend on both scattering length and 
lattice depth. Within first order perturbation theory, 
states of equal total quantum number Si + S2 + S3 + 11 
are coupled. To find the first-order energy corrections, 
we have to diagonalizc the matrix whose elements are 
given by 

<V4°> ff o + H'\4,f s ,) = 5 SS ,S SS ,E^ + <4°>'i4°l> 

(17) 

Here Sw is the Kronecker delta. Without interaction 
(a — > 0~) these states are degenerate as well, while for 
a being small and positive the interaction induced cou- 
pling becomes negligible, as was the case in the previous 
section. 

We exploit the symmetry of the perturbation term in 
order to analyze which relative and CM states will be 
coupled. This can be done conveniently for the case 
a = 0. Since the spatial symmetry of the wavefunction 
is not changed by the interaction, the same states will 
be coupled for a ^ 0, however, with coupling constants 
depending on a. The perturbation term is symmetric 
in cartesian coordinates, hence coupling occurs only be- 
tween states that have equal symmetry along the lattice 
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FIG. 3: (Color online) Total energy of the combined relative 
and CM system, as a function of the reciprocal scattering 
length. The upper thick lines represent the six states that 
are coupled by anharmonic terms. Thin gray lines denote the 
zeroth order HO levels. The scattering length is scaled on 
aho, with aho = i/2ft/ (mjjj) the harmonic oscillator length of 
the relative motion. 

axes. For a — > 0~ and without perturbation terms, the 
first excited energy consists of several degenerate states 
|V'ss) which can be divided into two sets. The first set 
consists of a CM ground state and excited relative states 
with labels 5 = {0,0,0} and s = {3,l,m}, and the sec- 
ond set consist of excited CM states and a relative ground 
state with S = {2,0,0} (and cyclic permutations), and 
,s = {0,0,0}. In the first case, quantum numbers can be 
I, m = {0, 0}, {2, 0}, {2, -2}, {2, +2}, hence four relative 
states in total. As required, the number of relative states 
that are coupled reduces to three when we take proper 
linear combinations of d-waves. So dictated by the sym- 
metry of the anharmonic terms, six states are coupled 
in total. Note that it is quite remarkable that a cou- 
pling to relative d-wave motion is possible, without an 
interatomic coupling on short range. This coupling be- 
tween s- waves and d-waves is indirect. The s- waves are 
coupled to excited CM states which in turn are coupled 
to d-waves, illustrating the long-range character of the 
anharmonic coupling. 

The main result of this section is shown in Fig. [3j To- 
tal energy for the combined relative and CM motion is 
shown as a function of the scattering length. Thick lines 
represent the six states that are coupled by anharmonic 
terms. Thin gray lines indicate the corresponding zeroth 
order HO levels. Also the perturbed uncoupled ground 
state is shown, using a thick solid line. The calculation 
is done for a lattice depth of Vq = 25E Ioc . Fig. [3] shows 
the presence of d-wave states in the upper two solutions, 
indicated with black dash-dotted and blue dashed lines. 
These solutions have only weak dependence on the s-wave 
scattering length, owing to the indirect coupling via the 
excited CM states. Below the d-wave states a dashed 
black solution is shown, which connects asymptotically to 
two consecutive levels, similar to the results of the Busch 
model for a s-wave scattering state. However, the three 
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FIG. 4: (Color online) Zoom in of Fig. [3] around resonance 
(1/a = 0). The scattering state (black dashed line) is cou- 
pled to one excited CM state (blue solid line), which con- 
tains a molecular bound state for positive a, while the other 
two states (dotted line) are degenerate and remain uncoupled 
around resonance. The splitting energy is indicated by -E ap iit- 

remaining solutions (dotted black, dotted blue and solid 
blue) are of a different nature. These reflect the presence 
of excited CM states, and asymptotically connect to the 
molecular state of the relative motion for positive a. A 
zoom-in around resonance (1/a = 0) is given in Fig. 2] 
It can be clearly seen that the coupling between relative 
and CM motion, caused by the long-range anharmonic 
terms, gives rise to a transfer of energy between these 
two motions. On resonance only the dashed and solid 
lines are coupled, giving rise to an energy splitting indi- 
cated by -Egpiit- The two dotted solutions are coupled at 
small negative values of the scattering length and with 
much smaller splitting energies. 

It is interesting to consider possible applications of this 
anharmonic coupling at long range, and interatomic cou- 
pling at short range. This can be done by exploiting 
different timescales when changing the scattering length. 
The scattering length can be changed by utilizing the 
magnetic field dependence of the scattering length via a 
Feshbach resonance. One can transfer for instance atoms 
from the lowest band into the next band by ramping 
the magnetic field, and by slowly ramping back associate 
molecules with an excited CM motion. This would re- 
sult in molecular energy levels that deviate significantly 
from the energy of ground state molecules labeled with 
5 = {0,0,0} and s = {1,0,0}. Such higher molecular 
levels correspond to molecules in different (partly-filled) 



Brillouin zones, which should be possible to detect [1 31 ] - 
We note that these excited CM molecules, compared to 
ground state molecules, have a larger spatial extent in 
the CM motion. 

The typical timescale r for slowly ramping back is 
given by r > h/E sp nt. In contrast to other lattice 
induced molecules, see e.g. [HI, these excited CM 
molecules do not arise from tunneling to neighboring 
sites, but from the anharmonic shape of a single lattice 
site. In addition, these excited CM molecules could be 
observable even for rather deep lattices (Vb > 50£r), 
since the anharmonic effects decay only weakly with in- 
creasing lattice depth, whereas tunneling effects decay 
exponentially. 



IV. CONCLUSION 

In conclusion, we proposed an easy to use method to 
solve for interacting atoms in optical lattices, where the 
relative and center-of-mass motion of the two interacting 
atoms arc coupled via the anharmonic terms of the lat- 
tice. The interactions are treated exactly using a bound- 
ary condition rising from a pseudo potential. The an- 
harmonic terms of the lattice potential are treated as a 
perturbation of the exact solution for two cold interact- 
ing atoms in a harmonic trap. We applied this model 
to the Mainz spin dynamics experiment 0, H| for / = 1 . 
The interaction energy is computed as the difference be- 
tween two-atom energy levels with and without interac- 
tion. This model gives a more rigorous interpretation 
of the experiment compared to previous descriptions, in 
terms of two-body scattering properties. We find that 
the derived scattering lengths agree within the experi- 
mental and theoretical error bars. Apart from applying 
our model to spin dynamics, we are in a good position 
to analyze future optical lattice experiments where the 
interactions are made very strong by utilizing Feshbach 
resonances. Strong interactions can induce coupling be- 
tween the relative and center-of-mass motion, which al- 
lows for an energy exchange between these two motions, 
and which can be used to produce (Feshbach) molecules 
with an excited center-of-mass motion. This model can 
also be used as a starting point for a description of pho- 
toassociation in an optical lattice near a Feshbach reso- 
nance [TBI. [l6j. 

This work was supported by the Netherlands Organi- 
zation for Scientific Research (NWO). 
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